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Abstract: We consider the semilinear wave equation in the radial case with conformal subcrit- 
ical power nonhnearity. If we consider a blow-up point different from the origin, then we exhibit a 
new Lyapunov functional which is a perturbation of the one dimensional case and extend all our 
previous results known in the one-dimensional case. In particular, we show that the blow-up set 
near non-zero non-characteristic points is of class C^, and that the set of characteristic points is 
made of concentric spheres in finite number in {-^ < < i?} for any R> 1. 
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1 Introduction 

We consider radial solutions of the following semilinear wave equation: 

(1.1) 



C/(0) = Uq and C/t(0) = t/i, 

where U{t) : a; G ^ C/(x, t) G M, C/q G Ri^^ ^ and Ui G Lf^^^^. 
The space Lj^^^ ^ is the set of all v in L^^^ such that 

\ 1/2 

v{x)\'^dx < +00, 




and the space HJ^^^^ = {v\v,\/v € Lf^c,u}- 
We assume in addition that 

4 

1< p and p < IH if iV > 2. (1.2) 

^ ^ - AT _ 1 - ^ ^ 

Since U is radial, we introduce 

u{r,t) = U{x,t) if r = \x\ (1.3) 
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and rewrite (1.1) as 

< dru{Q,t) = (}, (1.4) 
^ n(r, 0) = no(r) and ut{r,0) = ui{r), 

where u{t) : r £ — )■ u{r,t) £ M. We solve equation (1.1) locally in time in the space 
-'^loc ^ -'^k3c(^^) (^^^ Ginibre, Soffer and Vclo [6], Lindblad and Soggc [12]). Existence of 
blow-up solutions follows from ODE techniques or the energy-based blow-up criterion of 
Levine [11]. More blow-up results can be found in Caffarelli and Friedman [5], [4], Alinhac 
[1], [2] and Kichenassamy and Liftman [9], [10]. 

If u is a blow-up solution of (1.4), we define (see for example Alinhac [1]) a 1-Lipschitz 
curve r = {{r, T{r))} where r > such that the maximal influence domain D oi u (or the 
domain of definition of u) is written as 

D = {{r,t)\t<T{r)}. (1.5) 

r is called the blow-up graph of u. A point ro > is a non-characteristic point if there 
are 

So £ (0, 1) and to < T{ro) such that u is defined on Cro^T{ro),So ^{t> to} n {r > 0} (1.6) 

where C^, ^-j = {(r, t) \ t < t - S\r - r\}. We denote by 7^ C M+ (resp. S C M+) the set of 
non-characteristic (resp. characteristic) points. 

In a series of papers [16], [17], [20] and [18] (see also the note [19]), we gave a full 
picture of the blow-up for solutions of (1.1) in one space dimension. In this paper, we aim 
at extending all those results to higher dimensions in the radial case, outside the origin. 

Throughout this paper, we consider U{x, t) a radial blow-up solution of equation (1.1), 
and use the notation n(r, t) introduced in (1.3). We proceed in 3 sections: 

- in Section 2, we give a new Lyapunov functional for equation (1.4) and bound the solution 
in the energy space. 

- in Section 3, we study 7^, in particular the blow-up behavior of the solution and the 
regularity of the blow-up set there. 

- in Section 4, wc focus on S, both from the point of view of the blow-up behavior and 
the regularity of the blow-up set. 

2 A new Lyapunov functional and a new blow-up criterion 
in the radial case 

The extension of these results to higher dimensions in the radial setting involves the very 
beginning of our work in one dimension, namely the existence of a Lyapunov functional 
and the boundedness of the solution in similarity variables, performed in Antonini and 
Merle [3] and Merle and Zaag [13]. 

In [8] and [7] , Hamza and Zaag considered the following class of perturbed wave equa- 
tions for 

uu = /^u+\u\P-^u + f{u) + g{ut), (a;,t)GM^xM; (2.1) 
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where p satisfies (1.2) and for some q G [0,p), 

|/(x)| < C(l + and \g{x)\ < C(l + \x\). 

They showed in [8] and [7] that the argument of Antonini, Merle and Zaag in [3], [13], 
[15] and [14] extends through a perturbation method to equation (2.1). The key idea is 
to modify the Lyapunov functional of [3] with exponentially small terms and define a new 
functional which is in the same time decreasing in time and gives a blow-up criterion. 

It happens that the perturbation argument developed for equation (2.1) in [8] and [7] 
works for equation (1.4) with no further modification, as far as blow-up points different 
from the origin are considered. For the reader's convenience, we briefly recall the argument 
of Hamza and Zaag in the following. 

Given ro > 0, we recall the following similarity variables' transformation 

Wro{y,s) = {Tiro)-t) — u{r,t), y = s = -log(r(ro)-t). (2.2) 

The function w = satisfies the following equation for all y G (—1,1) and s > 
max (- logr(ro), - log ro): 

5> = ji:w- '^}^~^]l w + \w\P-'^w - ^^dsw - 2ydl,w + e-' ~ dyW, (2.3) 
(p — 1)^ p — I ro + ye~^ 

where jCw = -dy [p{l — y^)dyw) said p{y) = {1 — y^)p-^ . (2.4) 
Let us recall the Lyapunov functional in one space dimension 

which is defined in the space 

H = |g G Hi,, X LL(-1, 1) I ||g||^ = j'^ (g? + {q'.f (1 - y^) + pdy < +oo| . 

(2.6) 

Introducing 

F{w, s) = E{w) - e'" J wdswpdy, (2.7) 

we claim the following: 

Proposition 2.1 (A new functional for equation (2.3)) 

(i) There exists ^{p) > and So{N,p) G M such that for all ro > and for all s > 
max(-logr(ro),-S'o,-S'o - 21ogro), 

^^F{wr,{s),s) < j{p)e-'F{wr,{s),s) --^j^{dsWr,{s)f,^^dy. (2.8) 

(ii) (A blow-up criterion) There exists S\{p) G M such that ifW is a solution of equation 
(2.3) with ||I^(s)||/^p+i(B) locally hounded, and F{W{sq),so) < Q for some sq > Si{p), then 
W cannot be defined on the whole interval [so,oo). 
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Remark: From (i), we see that the Lyapunov functional for equation (2.3) is in fact 
H(wrQ{s), s) where 

H{w, s) = F{w, s)e-^(f , (2.9) 

not F{wrQ{s), s) nor E{'Wrf^{s)). 

Remark: Wc ahcady know from [13] and [15] that even in the non-radial setting, equa- 
tion (1.4) has a Lyapunov functional given by a natural extension to higher dimensions of 
E{wrQ{s)) (2.5). Unfortunately, due to the lack of information on stationary solutions in 
similarity variables in dimensions N >2, we could not go further in our analysis, and we 
had to stop at the step of bounding the solution in similarity variables. On the contrary, 
when = 1, we could obtain a very precise characterization of blow-up. 
Here, restricting ourselves to radial solutions, we find a different Lyapunov functional in 
higher dimensions (which exists even for supercritical p). Considering blow-up points dif- 
ferent from the origin, the characterization of stationary solutions in one space dimension 
is enough, and we are able to go in our analysis as far as in the one-dimensional case. 

Following our analysis in [13] and [16], we derive with no difficulty the following: 

Proposition 2.2 (Boundedness of the solutions of equation (2.3) in the energy 

space) For all ro > 0, there is a C2(ro) > and -S'2(ro) G M such that for all r G [^, ^] 
and s > S2{ro), 

£^ {{dyWr{s)f{l - 2/2) + (^,(,))2 + {dsWr{s)f + \Wr{sr+') pdy < C^ivo). 

Proof: The adaptation is straightforward from our analysis in [13] and Proposition 3.5 page 
66 in [16]. The only difference is in the justification of the limit at infinity of E{wr^-,{s)), 
which follows from the limit of H(wro{s), s) defined in (2.9). In fact, we know from 
Proposition 2.1 that H{wrQ{s),s) is decreasing and bounded from below, and such an 
information is unavailable for E{wro{s)). 

Proof of Proposition 2.1: 
(i) Consider ro > 0, s > max(— logr(ro), — log ^) and write w = for simplicity. From 
the similarity variables' transformation (2.2), we see that 



ro + ye G 



ro 3ro 
2 ' 2 



(2.10) 



Multiplying equation (2.3) by dgwp and integrating for y G (—1,1), we see by definition 
(2.5) of E{w) that 



^E{w{s)) = i-- / {dMs)fT-^dy + {N-l)e-^ [ dMs)dyw{s)% (2.11) 

where r is defined in (2.10). Using (2.10), we write 



/I 
^dsW{s)dyW{s)% 



(2.12) 

< ( ^ ) I _^idyw{s)rpil - y-)dy+ ^-f^ J [^—-^ j JJdsw{s)r^^dy. 



'p-^\ ^ ^^2 n , ^ 4 \fN-iyf\^ ,,,2 p 
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Recalling the following Hardy-Sobolev estimate (see Appendix B page 1163 in [13] for the 
proof): 

/ h^-^dy<C f h^pdy + C [ {h'{y)fp{l-y^)dy, (2.13) 

7-1 ^ — y J-i j-i 

we use the notation I{s) = J dswwpdy and write from equation (2.3) for any e > 
and s > - log (2(fe) ) > 

e'l'is) > j'^\w{s)\P+'pdy-{l + e)j'jdyw{s)\^p{l-y^)dy (2.14) 



e 



■ \dsw{sr^dy - (e + ^^) Jjw{s)fpdy. 



Using (2.11), (2.12) and (2.14), we get that (2.8) follows by definition (2.7) of F{w,s), 
provided that we fix e = e{p) > small enough and take s large enough. This yields (i) of 
Proposition 2.1. 

(ii) If is a solution of equation (2.3), then we write by definition (2.9) of H: 

-7(^)6-" rl 



HiW{s),s) > — — / \W{s)\P+'pdy 

J-i 



+e-7we- ( { - - — J jjdsWis)rpdy+^j^—^ ~^)J p'^y 



^ I ms)r'pdy 

if s > Si{p) for some Si{p) G M large enough. Using this inequality together with the fact 
that H(W(s),s) is decreasing by the remark following Proposition 2.1, we see that the 
argument used by Antonini and Merle in Theorem 2 page 1147 in [3] for the equation (1.1) 
works here and we get the blow-up criterion. This concludes the proof of Proposition 2.1. 



3 Blow-up results related to non-characteristic points 

Let us first introduce for all |d| < 1 the following solitons defined by 

('1_^2\^ /2(p + l)\p^ 
K{d,y) = Ko ^ where Ko = 7 ttt and |y| < 1. (3.1) 



(l + dy)^ \iP-^ 

Note that K{d) is a stationary solution of (2.3) in one space dimension. 
Adapting the analysis of [16] and [17], we claim the following: 

Theorem 1 (Blow-up behavior and regularity of the blow-up set on TZ) 
(i) (Regularity related to TV) TZ ^ ^ , Tlf] Wj^ is an open set, and x ^ T{x) is of class 
on 7^n Ml. 
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(ii) (Blow-up behavior in similarity variables) There exist /xq > and Cq > such 
that for all ro £ TZCi U*^, there exist 6{ro) = ±1 and so('^o) > — log T(ro) such that for all 
s > So- 

< Coe'-i'°^'-'°\ (3.2) 

H 

Moreover, E{wro{s)) — )• E(ko) as s ^ oo. 

Remark: If € 7^, the asymptotic behavior of wq remains open. 

Proof: The proof is the same as for the one-dimensional case treated in [16] and [17], 
with some minor adaptations. For that reason, we don't give the full proof here. We will 
instead ask the reader to follow our proof in those papers, and provide him only with the 
delicate points he may face in adapting the proof to the higher-dimensional radial case. 
As in the one-dimensional case, we proceed in two steps: we first prove (ii) with a param- 
eter do(''o) instead of T'(ro), then we prove (i) and the fact that in (ii) do{ro) = T'{ro). 

Step 1: Proof of (ii) with do(ro) instead of r'(ro) in (3.2) 

This statement is the twin of Corollary 4 page 49 given in [16] in the one-dimensional 
case. 

Consider ro G 7?. fl W^. As in that case, the proof has two major steps: 

- Approaching the set of stationary solutions, in the sense that for some 9{ro) = ±1 
and do{ro) G (0, 1), we have 

inf \\wro{-,s) - e{ro)K{d,-)\\HU-i,i) + \\dsWj;jL2(^_iu ^ (3.3) 

\d\<do{ro) 

as s — )• OO (see Theorem 2 page 47 of [16] for the statement in one space dimension). 
Note that such a statement is still unavailable for non radial solutions of equation (1.1), 
since we have no classification for stationary solutions of the multi-dimensional version of 
equation (2.3). Fortunately, in the radial case, we naturally see from equation (2.3) and 
the estimate (2.10) that we only need the classification in one space dimension, given in 
Proposition 1 page 46 in [16]. This is the reason why we restrict ourselves to the radial 
case in this paper. As for the proof of (3.3), the reader has to see Section 3 page 60 in 
[16]. The only delicate point is the adaptation of the proofs of Lemma 3.3 page 62 and 
Lemma 3.4 page 64. Indeed: 

- In the proof of Lemma 3.3, we need a Duhamel formulation for the one-dimensional 

semilinear wave equation (see estimate (63) page 63 in [16]). This formulation has to be 
replaced by the radial version, which may be derived from Shatah and Struwe [21]. 

- As for the proof of Lemma 3.4, due to the fact that equation (1.4) is no longer invariant 
under the scaling 

A ^ uxi^, t) = X^u{X^, Ar), 

we need to understand the continuous dependence of the solutions of the following family 
of equations 

d?u = dju + ^^~]^^ dru + \u\^-^u, (3.4) 

with respect to initial data and the parameters a; > and A > (including the limit as 
A ^ 0), and this is a classical estimate. Apart from this point, there is no other problem 
in the adaptation. 



6 



- Proof of the convergence to a 0{ro)K{do{ro)) for some (io(?'o) € (—1,1). From (3.3) 
and the monotonicity of the Lyapunov function H{wrQis), s) (2.9), we see that (h) (with 
a parameter do(ro) instead of T'{ro) in (3.2)) follows directly from the following trapping 
result. Note that the statement of this result is different from the analogous result in one 
space dimension given in Theorem 3 page 48 in [16]. This difference is due to the fact that 
the equation (2.3) in similarity variables depends on a parameter tq > and contains a 
new term of order (it is no longer autonomous). For a further purpose, we give in the 
following the radial case's version of the trapping result, valid uniformly for all ro > po 
for some po > 0: 

Theorem 2 (Trapping near the set of non zero stationary solutions of (2.3)) 

For all pq > 0, there exist positive eg, po and Cq such that for all e* < eg, there exists 
so(e*) such that if ro > po, s* > sq and w G C{[s* , (x>),T-l) is a solution of equation (2.3) 
with 

ys>s*, E{w{s))>E{Ko)-e-i, (3.5) 

and 

<s*) \_,.( <d*,-) 

dMs*) J V 



< e* 



H 



for some d* G (—1,1) and a;* = ±1, then there exists d^o G (—1,1) such that 



largthdoo — argthd*| < Coe*, 



and for all s> s* , 



w{s) 
dsw{s) 



K{do 





n 



Remark: The function argth is the inverse of the hyperbolic tangent function. It is given 
by argth d = ^ log (^i^s) • Theorem 2 holds under a weaker condition where we replace 
the e~2 in (3.5) by C*e~^ for some constant C*(po) > 0. 

Proof. The proof can be adapted from the proof of Theorem 3 in [16] given in Section 5 
page 99 in that paper. Up to replacing u by —u, we may assume u* = 1. As in [16], we 
linearize equation (2.3) by introducing 



q{y,s) = w{y,s) - K{d{s),y) 
where the parameter d{s) G (—1, 1) is chosen by modulation so that 



(3.6) 



(3.7) 



and ttq is the projector on the null mode of the linearized operator. Then, we decompose 
q into two parts, according to the spectrum of the linearized operator as follows: 

- its projection on the mode A = 1, whose norm is bounded by some |ai(s)|; 

- its projection on the negative part of the spectrum, whose norm is bounded by some 
a_(s). 

Thanks to this decomposition, we write 



1 

Co 



(|ai(s)| + a_(s)) < \\qis)\\n<Coi\aiis)\+a.is)) 
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for some Cq > 0. 

Projecting the linearized equation according to this decomposition, we see that we have 
exponentially small perturbations coming from the new term in the similarity variables 
equation (2.3), which is uniformly bounded since the parameter ro satisfies rg > po > 0, 
in the sense that 



(iV-l)e-% 

o-,,w 



ro + ye 



-s 



< —{N - l)e-'\dyw\ as soon as s > - log^. (3.8) 
Po 2 



More precisely, we have the following statement, which is the new version of Proposition 
5.2 page 103 in [16]: 

Proposition 3.1 There exists £2 > such that ifw a solution to equation (2.3) satisfying 
(3.7) and 

Ms)\\h<^ 

at some time s for some e < 62, where q is defined in (3.6), then: 
(i) (Control of the modulation parameter) 

\d'\ < Co(l - d'^){af + a J) + Ce-'{1 - d^). 
(a) (Projection of the linearized equation on the different eigenspaces of L^) 



i' 

ai — a\ 



< Co [a\ +a-^)+ Ce 



-s 



for some R-{s) satisfying 

\R-{s)\ < Co{a\^ + a-^)^ where p = mm{p, 2) > 1. 
(Hi) (Additional relation) 



4:1 QiQ2p<-iaJ^ + Co I g:i,2T^ + Co^r + Ce- 



;2 



ds J_i 5 J_i ' 1 - y 

(iv) (Energy barrier) If moreover (3.5) holds, then 

a\{s) < Coa-{s) + Ce-'/^. 

Proof. In comparison with the one dimensional case, the linearized equation has an ad- 
ditional exponentially small term. Therefore, the adaptation consists in projecting that 
term on the different components, and there is no difficulty for this. The projection of the 
remaining terms is the same as the one-dimensional case. See the "Proof of Proposition 
5.2" page 104 in [16] for the one-dimensional case. This concludes the instructions for the 
adaptation of Proposition 3.1. ■ 

Now, with this proposition, there is no particular difficulty in deriving Theorem 2 as in 
the one-dimensional case. See Section 5.3 page 113 in [16]. This concludes the instructions 
for the adaptation of Theorem 2. ■ 
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Step 2: Proof of (i) and the fact that (ii) holds with a parameter do{ro) = 
r'(ro) in (3.2) 

Note that the fact that 7?. 7^ follows by the same argument given in the remark 
following Theorem 1 page 58 in [17]. 

The corresponding statement to this step in the one-dimensional case is Theorem 1 page 58 
in [17]. As in that case, we need a rigidity Theorem (or a Liouville Theorem) in similarity 
variables. Let us insist on the fact that we don't need any adaptation for the Liouville 
Theorem and that the version we need is indeed the one-dimensional version stated in 
Theorems 2 and 2' in pages 58 and 59 of [17]. While adapting the proofs to the radial 
case, the reader should pay attention to two facts: 

- Again, we need to use the continuity of solutions of equation (3.4) with respect to initial 

data and to the parameters (x. A). 

- The functional E{w) (2.5) is no longer nondecreasing. Fortunately, we can replace it 
by the Lyapunov functional H(w, s) (2.9) which is decreasing. The error we make in this 
replacement is exponentially small as one sees from (2.7), (2.9) and Proposition 2.2. 

- We need the trapping result here too. Of course, the new version given in Theorem 2 
applies. 

Apart from these three remarks, the adaptation is straightforward. One has just to 
follow the proof of Theorem 1 given in Section 2 page 60 in [17]. 

4 Blow-up results related to characteristic points 

In Proposition 1 in [20], we showed the existence of a solution to equation (1.1) in one space 
dimension such that 5 7^ 0. Artificially adding other coordinates, this one-dimensional 
solution can be considered as a multi-dimensional solution with <S 7^ 0. Adapting the 
argument of [20] to the radial case, we prove the following: 

Proposition 3 (Existence of radial solutions with a non zero characteristic 
point) There exists {uo,Ui) such that the corresponding solution of equation (1.4) has 
a non zero characteristic point. 

Remark: In this case, the multi-dimensional version U{x,t) = u{\x\,t) has a sphere of 

characteristic points. 

In fact. Proposition 3 is a consequence of the following: 

Theorem 4 (Existence and generic stability of characteristic points) 

(i) (Existence) Let < ai < 02 be two non characteristic points such that 

^oi(s) — >■ 0{ai)K{dai, •) as s 00 with 9{ai)6(a2) = —1 

for some da^ in (—1,1), in the sense (3.2). Then, there exists a characteristic point 
c £ (01,02)- 

(ii) (Stability) There exists eo > such that if \\{Uo, Ui) — (C/q, J7i)||hi x (r^) ^ ^o? 

loc.ii loc.u ^ ' 

then, u{r,t) the solution of equation (1.4) with initial data {uo,ui)(r) = (Uo,Ui){x) if 
r = \x\ blows up and has a characteristic point c G [01,02]. 
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Remark: This statement is different from the original one (Theorem 2 in [20]) by two 
natural small facts: we take positive points ai and 02 in (i), and we use the multi- 
dimensional norm in (ii) (of course, from the finite speed of propagation, it is enough 
to take a localized norm instead). 

Indeed, let us first derive Proposition 3 from this theorem, then prove this latter. 

Proof of Proposition 3 assuming Theorem 4- From the finite speed of propagation and 
the uniqueness of the solution to the Cauchy problem, one can take {uo,ui) with large 
plateaus of opposite signs so that for some < ai < 02 and large T > 0, the solution 
remains space independent in the backward cones {r — < T — t}, conserves its sign 
there, and blows some time T{ai) < T. Since the solution is space independent 

around Oj, the blow-up time is locally constant and the point Oj G 7^ with T'{ai) = 0. 
Using the description of the blow-up behavior in the non-characteristic case stated in (ii) 
of Theorem 1, we see that (y, ,s) — > 9{ai)K as s — > cx) with 6(01)9(02) = — 1, since the 
plateaus have opposite signs. Therefore, the hypotheses of (i) of Theorem 4 are fulfilled 
and we have the desired conclusion of Proposition 3. ■ 

Now, we give indications on the adaptation of the proof of Theorem 4 from the one- 
dimensional case. 

Proof of Theorem 4- There is no difficulty in adapting to the present context the proof 
of Theorem 2 of [20] given in Section 2 of that paper, except may be for the continuity 
of the blow-up time with respect to initial data, stated in Proposition 2.1 of [20], where 
some natural extensions to the radial case are needed. ■ 

We also have the following result which relates the existence of characteristic points 
to the sign-change of the solution: 

Theorem 5 (Non-existence of characteristic points if the sign is constant) Con- 
sider u(r,t) a blow-up solution of (1.4) such that u(r,t) > for all r G (ao,bo) and 
to <t < T(r) for some real < oq < 60 o,i^d to > 0. Then, (ao, bo) C TZ. 

Remark: This statement is exactly the same as the original (Theorem 4 in [20]). In 

particular, it is valid with oq = 0. 

Proof. This result follows from Theorem 6 below exactly as in one space dimension. See 
the proof of Theorem 4 given in Section 4.1 in [20]. ■ 

Now, given tq G 5 H M*^, we have the same description for the asymptotics of tt;^,, as 
in the one-dimensional case. More precisely, the following holds in the radial case, outside 
the origin (for the statement in one space dimension, see Theorem 6 in [20]): 

Theorem 6 (Description of the behavior of where ro is characteristic) Con- 
sider u(r,t) a blow-up solution of (1.4) and ro G 5 n M^. Then, it holds that 

and E(wro(s)) k(ro)E(Ko) (4.1) 

H 
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as s ^ oo, for some 



k{ro) > 2, (4.2) 
e* = eU-iy+' (4.3) 

and continuous di{s) = — tanh(^j(s) G (—1,1) for i = 1, k{ro). Moreover, for some 
Co > 0, for all i = 1, k{ro) and s large enough, we have 

0„_(,_(Mrci±i))(P^,„,,|<C„. (4.4) 

Proof. As in the one-dimensional case, the proof of the asymptotic behavior and the 
geometric results on S (see Theorem 8 below) go side by side. For that reason, we leave 
the proof after the statement of Theorem 8. ■ 

Extending the definition of /c(ro) defined on S in Theorem 6 by setting 

Vro G 7^, k{ro) = 1, 

we get the following result on the energy behavior from the asymptotic behavior at a 
non-characteristic point (see (ii) of Theorem 1) and at a characteristic point (see Theorem 
6): 

Corollary 7 (A criterion for non characteristic points) 

For all ro > 0, there exist Cs{ro) > and Ss{ro) G M such that: 

(i) For all r G [^, ^] and s > S^, we have 

Eiwris)) > kir)E{Ko) - C3(ro)e-^ 

(ii) If for some r G ^] and s > S3, we have 

Eiwris)) <2E{Ko)-C3{ro)e-', 

then r eTZ. 

Remark: With respect to the one-space dimension statement (Corollary 7 in [20]), this 
statement has additional exponentially small terms. This comes from the fact that the 
functional E{w) is no longer decreasing, and that one has to work instead with the func- 
tional H{w, s) (2.9) which is decreasing, and differs from E{w) by exponentially small 
terms, uniformly controlled for r G [^, ^] thanks to the uniform estimates of Proposi- 
tion 2.2. 

Proof: If one replaces E(w) by H{w, s), then the proof is straightforward from Theorems 1 
and 6 together with the monotonicity of H{w, s) (see (2.9) and (2.1)). Since the difference 
between the two functionals is exponentially small, uniformly for r G [^,^] (see (2.9), 
(2.8) and Proposition 2.2), we get the conclusion of Corollary 7. ■ 

Finally, we give in the following some geometric information related to characteristic 
points (for the statement in one space dimension, see Theorem 1, Theorem 2 and the 
following remark in [18]): 
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Theorem 8 (Geometric considerations on 5) Consider u{r,t) a blow-up solution of 
equation (1.4). 

(i) (Isolatedness of characteristic points) All characteristic points different from the 

origin are isolated. 

(ii) (Corner shape of the blow-up curve at characteristic points) If r^ G 5 n 

and < |r — ro| < then 

< T'{r) + ^ < 



Co|log(r -ro)| 2 I' '0| |log(r-ro)| 2 

for some > and Co > 0, where k{ro) > 2 is the integer defined in (4.2). 
Proof. See below. 

Remark: Integrating the estimate in (ii) of this theorem, we see that 

< T{x) - T{xo) + \x- xo\ < °' (J„U)(.-i) - (4-5) 

Col log(a; - xo)| 2 |log(x-xo)| 



2 



Remark: Note from (i) that the multi-dimensional version U{x,t) = u{\x\,t) has a finite 
number of concentric spheres of characteristic points in the set {-^ < |a;| < i?} for every 
R> 1. This is consistent with our conjecture in [18] where we guessed that in dimension 
N > 2, the (N — l)-dimcnsional Hausdorfi^ measure of S is bounded in compact sets of 
R^. Note that this conjecture is related to the result of Velazquez who proved in [22] 
that the {N — l)-dimensional Hausdorff measure of the blow-up set for the semilinear heat 
equation with subcritical power nonlinearity is bounded in compact sets of ' 



As a consequence of our analysis, particularly the lower bound on T(r) in (4.5), we 
have the following estimate on the blow-up speed in the backward light cone with vertex 
(ro,T(ro)) where ro > (for the statement in one space dimension, see Corollary 3 in 
[18]): 

Corollary 9 (Blow-up speed in the backward light cone) For all ro > 0, there 
exists C4(ro) > such that for all t € [0, T(ro)), we have 

|log(TM-t)|^^ sup M.,t)|<^^('-»'l'°«<^<''»'-*'l'^ 



2 ^ ^'^f |"'V^'''yi — 2_ 

-1 



C4(ro)(r(ro) - t)f-i |a;-ro|<r(ro)-t {T{ro) - t)p- 

Remark: Note that when r^ £ TZCi M!j_, the blow-up rate of u in the backward light cone 
with vertex (ro,r(ro)) is given by the solution of the associated ODE u" = vP. When 
ro G 5 n R;^, the blow-up rate is higher and quantified, according to A;(ro), the number of 

solitons appearing in the decomposition (4.1). 

Proof. When ro G TZ, the result follows from the fact that the convergence in (3.2) is true 
also in L°° x from (3.3) and the Sobolev embedding in one dimension. When ro G 5, 
see the proof of Corollary 3 of [18] given in Section 3.3 of that paper. ■ 

Proof of Theorems 6 and 8: The proof follows the pattern of the original proof, given 
in [16], [20] and [18]. In the following, we recall its different parts. 
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Part 1: Proof of (4.1) without (4.2) nor (4.3) and with the estimate 

Cj+i(s) - Ci(s) ^ oo as s ^ oo (4.6) 

instead of (4.4) (note that both (4.4) and (4.6) are meaningful only when fc(ro) > 2). 

The original statement of this part is given in Theorem 2 (B) page 47 in [16] and 
the proof in section 3.2 page 66 in that paper. Note that this part doesn't exclude the 
possibility of having k{ro) = or k{ro) = 1. The adaptation is straightforward. As in the 
non-characteristic case above, one has to use the Duhamel formulation in the radial which 
may be derived from [21]. 

Part 2: Proof of (4.3), (4.4) and the upper bound in (4.5), assuming that 
(4.2) is true. 

The original statement is given in Propositions 3.1 and 3.13 in [20]. The reader has to 
read Section 3 and Appendices B and C in that paper. The adaptation is straightforward, 
except for the effect of the new term in equation (2.3), which produces exponentially small 
terms in many parts of the proof In particular. Lemma 3.11 of [20] has to be changed by 
adding Ce"^ to the right of all the differential inequalities. 

Part 3: Proof of (4.2) and the fact that the interior of S is empty. 

The original statement is given in Proposition 4.1 of [20]. Here, the adaptation is 
not only delicate, but we need a new argument to rule out the occurrence of the case 
where, locally near the origin, the blow-up set of the multi-dimensional version U{x, t) is 
a forward light cone with vertex (0,T(0)) (we actually prove a stronger result, see Lemma 
4.5 below). For this reason, we will give the good version of Proposition 4.1 in [20] and 
outline its proof in the radial case. More precisely, we claim the following: 

Proposition 4.1 

(i) The interior of S is empty. 

(ii) For all ro € 5 n W^, k{ro) > 2. 

Remark: Please note that in (i), the information is about S, whereas in (ii), we have to 
restrict to 5nM+. 

As in [20], this proposition is a consequence of the following Lemmas, which we restate, 
since some statements surprisingly remain valid even at the origin, whereas others are valid 
only outside the origin: 

Lemma 4.2 (Characterization of the interior of 5) For any < ri < r2, the fol- 
lowing statements are equivalent: 

(a) (ri,r2) G S. 

(b) There exists r* G [ri, r2] such that for all r G [ri, r2], T{r) = T{r*) — |r — r*\. 

Remark and Proof. Note that this lemma is valid also at the origin. The proof is the same 
as for Lemma 4.2 in [20]. ■ 

Lemma 4.3 Consider < ri < r2 such that e = _ _|_-|^ Then, 

(i) for all r G [ri, r2], T{r) = T{ri) + e{r - ri), 

(ii) (ri,r2) G S. 
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Remark and Proof: Note that this lemma is valid at the origin too. The proof can be 
adapted straightforwardly from the proof of Lemma 4.3 in [20]. ■ 

Lemma 4.4 (Boundary properties of 5) 

(i) For all To e dS D R^, k{ro) 7^ 0. 

(ii) Consider E dS Ci with k{ro) = 1. If there exists a sequence r^j G 7^ converging 
from the left (resp. the right) to ro, then vq is left-non-characteristic (resp. right-non- 
characteristic). 

Remark: Unlike the two previous lemmas, this lemma is valid outside the origin. This 
is due to the fact that we strongly need the structure in similarity variables, which is 
available only outside the origin. We mean by ro is left-non-characteristic (resp. right- 
non-characteristic) that it satisfies condition (1.6) only for r < ro (resp. for r > ro). 
Proof. The only delicate point in the adaptation of the proof from the proof of Lemma 
4.4 in [20] is in Claim 4.5. Indeed, we need to choose there the time i close enough to 
r(a;o) so that we can apply the trapping result stated in Theorem 2. Remember that this 
restriction in the trapping result comes from the fact that the equation (2.3) is no longer 
autonomous, and that the new term in (2.3) becomes small when s is large (see (3.8)). 
We would like to add that Claim 4.5 of [20] and its proof given in Appendix D don't use 
the equation satisfied by so the proof is rigorously the same. ■ 

In addition to the above lemmas, we have to add a new ingredient to the proof: the 
blow-up set of the multi-dimensional version U{x, t) is always strictly under the forward 
light cone of vertex (0, r(0)). More precisely, we make the following statement: 

Lemma 4.5 (The blow-up set is strictly under the forward light cone with 
vertex (0,T(0))) For all r > 0, we have T{r) < r(0) + r. 

Let us first use the previous lemmas to derive Proposition 4.1, then we will prove 
Lemma 4.5. 

Proof of Proposition 4.I assuming that Lemma 4-5 is true: 
(i) Arguing by contradiction, we assume that S fl M!^ contains an open non empty interval. 
Maximizing that interval, we can assume that a maximal interval (a, 6) is included in 
5 n M+ with < a < 6 < +00. Prom Lemma 4.2 and the fact that for all r > 0, T{r) > 0, 
we have two cases: 

Case 1:6 = +00 and for all r > a, T(r) = T{a) + r — a; 

Case 2: < 6 < -hoc, b e dS and for all r G [a,b], T{r) = T{c*) - |r - c*\ for some 
c* G [a, 6]. 

Let us find a contradiction in these two cases, assuming first that a = then a > 0. 

- If a = 0, then we see from Lemma 4.5 that a contradiction follows in Case 1 or in 
Case 2 if c* > 0. Now, if Case 2 holds with c* = 0, then, 

Vr G [0, 6], T(r) = r(0) - r and 6 G n dS. (4.7) 

If k{h) = 0, then a contradiction follows from (i) of Lemma 4.4. 

li k{h) = 1, then from the fact that 6 G 95, there exists a sequence rn & Tl converging to 
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b. Since (0, b) £ S hy hypothesis, we have rn > b for n large enough. Using (ii) of Lemma 
4.4, we see that b is right non-characteristic. Since 6 is clearly left non-characteristic from 
(4.7), this means that b eTZCi R*^. Since TZ fl R*^ is open from (i) of Theorem 1, this is in 
contradiction with the fact that b e dS. 

If k(b) > 2, then we know from Part 2 above that the upper bound in (4.5) holds, which 
means that the blow-up set is corner shaped near b, and this is a contradiction by (4.7). 

- If a > 0, then, a G dS since the interval (a, b) is maximal. If Case 2 holds with 
c* = a, then the proof is exactly the same as in the case given above where a = and 
Case 4 holds with c* = 0. Now, if Case 1 holds or Case 2 with c* > a holds, then 



since the interval (a, b) is maximal. The situation is symmetric with the situation where 
a = and Case 2 holds with c* = 0, and where we found a contradiction at the point b. 
Here, the contradiction follows in the same way, but at the point a. This concludes the 
proof of (i) in Proposition 4.1. 

(ii) The proof is exactly the same as the proof of the analogous statement in [20] (see 
The proof of Proposition 4.1 in Section 4.1 of that paper). 

This concludes the proof of Proposition 4.1 assuming that Lemma 4.5 is true. It remains 
then to prove Lemma 4.5. ■ 

Proof of Lemma 4-5: Prom invariance by time translation, we may assume that 



We proceed by contradiction and assume that for some ro > 0, we have T{ro) > vq. 
Recalhng that r r(r) is 1-Lipschitz, we see that T{ro) = vq. Using Lemma 4.3 (which 
is valid when ri = 0), we see that 



Recalling that p is subcritical or critical with respect to the existence of the conformal 
invariance (see (1-2)), we have the following bound on the blow-up rate from Theorem 1 
page 1149 in [13] and Theorem 1 page 397 in [15]: 



Vr G [a, c*], r(r) = T(a) + r - a and a eR^n dS, 



r(o) = 0. 



Vr G (0, ro), r(r) = r and (0, ro) C S. 



(4.8) 



In particular, by definition (1.6) of a non-characteristic point, 



G 7^. 



Vs > So 



^ < \\Wo{s)\\h^^\y\<i) + ||9.Wo(s)|U2(|;^|<i) < K, 



(4.9) 



where so € M, -fC > and 



VyGM^, Wo{Y,s) = wo{\Y\,s). 



(4.10) 



We claim that in order to conclude, it is enough to prove that 




(4.11) 
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and eo > 0. Indeed, if (4.11) holds, then we write from (4.9), 

Vn G N, jy \woiy, Sn)\P+'dy > | (4.12) 

for some (5 G (0, ^). 

Following our argument for the proof of Claim 3.12 page 77 in [17], wc take b £ (0, rg). 
Recalling from (4.8) that T(b) = b, we write by definition of similarity variables (2.2) that 

My', s') = (1 - be^'r^^woiy, s) with y = ^ ^ and s = s' - log(l - be''). (4.13) 

1 — oe* 

Introducing 

4 = - log(6 + e-""), y[{s') = 5-b{l + 5)e'' and y^(s') = 1 - 5 - 6(2 - 5)e'' (4.14) 
we see that 

< e*"6 < 1 and e*"6 ^ 1 as n ^ oo, (4-15) 

hence 

j/^(4)G(-l,l), 2/i(s) = 5, y2{s) = l-5, yy' &{y'MMs'n)l p{y')>^—^^^. 
Therefore, we write from (4.13) and (4.12), 



/ \w,{y',s'^)r^p{y')dy'> / \w,{y' , s'^)^' p{y')dy' 

(l~be<)~^ f^-\ , , eo(l-6e<)"i^ 
> Fo(y, Sn)r dy > ^^^^ )• oo, as rn- GO 

from (4.15). This contradicts the bound 

ljw,{y',s'jr'p{y')dy'<Co{b) 

stated in Proposition 2.2. Therefore, it is enough to prove (4.11) in order to conclude the 
proof of Lemma 4.5. 

Let us proceed by contradiction in order to prove (4.11), and assume that 

f \Wo{Y, s)\P^^dY ^ as s ^ oo. (4.16) 
7|y|<i 

Therefore, we see from (4.9) that for s large enough, we have 
1 

Introducing for all n G N, 



77^ < ||VW'o(5)||L^(|y|<i) + \\dsWois)\\L2^\Y\<i) < K. (4.17) 



Vnii, r) = (1 - t)-^^Wo{Y, s) with Y = s = n- log(l - r), (4.18) 

1 — r 
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we see from the definitions (4.10), (2.2) and (1.3) of Wq, wq and u that Vn is radial in the 
sense that 

Vn{^,T)=Vnm,T) (4.19) 

and that Vn is a solution of the multi-dimensional equation (1.1) in the backward fight 
cone of vertex (0, 1) above the section at time r = 0, in the sense that 

VrG[0,l), V|e|<l-r, d^Vn = AVn + \Vnr^Vn. 

Using (4.18), (4.17) and (4.16), we see that for n large enough and for all r G [0, 1), 



jV 2(p+l) 



p-1 



JV-2(P+1) 

^'^^LsJ^ < ll^-^nWIIi^dCKi-.) + l|VK(r)||i.(|^|<,_,) < 16KH1 - r) 

N i 

l|Ki(T)||LP+i(|^|<i-r) + llK(r)||L2(|^l<i_^) < e„(l-r)p+i f-i 

(4.20) 

where e„ ^ as n ^ oo. Since N — ^^^^p- < —2 from the condition (1.2) on p, there 
exists 

To{K)e {0,1) 



such that 

l|5.K(ro)||i.(|^l<i_,„) + ||VK(To)||i2(|^|<i_,„) > 
Applying the following one-dimensional Sobolev inequality 

/(l)'<C||/||^.(i,i)||/||^,(i,i) 
to the radial version Vn (4.19), we see from (4.21) and (4.20) that 

sup Vn{C,Tf<C{K) sup ||K||l2(|^|<i_^)||14||^i(|^I<i_^) < C(iC)e„ 

(e,r)6Bro 0<r<ro 



(4.21) 

(4.22) 



(4.23) 
(4.24) 



where 

^ro = {(e,r)|0<T<To, |ei = l-r} 

is the lateral boundary of the portion of the backward light cone of vertex (0, 1) located 
between the sections at r = and t = tq. Introducing the following local energy defined 
in the section of the backward light cone of vertex (0, 1) 



J\e\<i-7 



'iei<i-T L 

we write from classical estimates. 



{drVn{^,r)f + (VF„(C,r))2 - 



p+1 



S{Vn{To)) - S{Vn{0)) 



'V2 



B. 



p+i 



da 



V2(p+i: 



-da 



where B^^ is defined in (4.24). Using (4.20), (4.21) and (4.23), we see that for n large 
enough, we have 



l|5rK(ro)||i2(|^l<i_,„) + ||VK„(ro)||i2(|^,<i_,„) 

2(|^|<l) + l|V^n(0)llL2(|^|<l) 



< \\drVnm\hm<l) + l|VF„(0)||22(,„^,, + < 17i^2 
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This is a contradiction by (4.22). Thus, (4.11) holds and Lemma 4.5 is proved. Since 
Proposition 4.1 follows from Lemma 4.5, this concludes the proof of Lemma 4.5 and 
Proposition 4.1 too. ■ 



Part 4: Proof of Theorem 8 

The analogous statement in one space dimension is given in Theorems 1 and 2 in [18]. 
Thus, we need to say how to adapt the analysis of the paper [18] to the radial case. Let us 
recall the strategy of the proof from Section 1.3 in that paper. Consider u{r,t) a blow-up 
solution of equation (1.4) and vq £ S (1 The decomposition of ^^^(y, s) is given in 
Theorem 6 (up to replacing u{r,t) by —u{r,t), we may assume that = —1). 

To prove that tq is an isolated characteristic point, the only tools we have are the 
energy criterion in (ii) of Corollary 7 and the trapping result stated in Theorem 2. Note 

that due to the fact that equation (2.3) is no longer autonomous and depends on the 
considered blow-up point, the uniform version of the trapping we stated in this paper is 
strongly needed (we take po = ^ in Theorem 2). 

In order to use these tools, we have to find the behavior of Wr for r near tq. A simple 
idea is to start from the decomposition (4.1) for and the fact that the blow-up set is 
locally different from a straight line (it is in fact corner shaped by the Tipper bound in 
(4.5) already proved in Parts 3 and 2 above), and use the transformation (2.2) first to 
recover the behavior of u{r,t), then the behavior of Wr{y, s) for r near vq. Two problems 
arise in this simple idea: 

- we can't have information on Wr{y,s) for all y G (—1,1), since information on the 
whole interval (—1,1) would involve information on Wr(^{y,,s) for \y\ > 1, and this is 
unavailable (at least at time s) because of the finite speed of propagation; 

- the relation between and Wr we get from (2.2) depends explicitly on the value 
of r(r) which is an unknown. The value of T{r) specified by the upper bound in (4.5) 
changes the range of s for which we have information. 

To overcome these problems, we proceed in 3 steps: 

Step 1: Initialization of the behavior of Wr{s) 

Here, we use (4.1) and continuity arguments to show that for r close enough to vq and 
s = Lfn-\-i large enough, Wr is close to a sum of k solitons 

k 

^i-iyKl{di{Lm+l),i>iiLm+l)) 

1=1 

where 

Klid,i^,y) = no^ (4.25) 

{l + dy + iy)p-i 

Note that K^(d, ±e*) are heteroclinic orbits of the one-dimensional version of equation 
(2.3) connecting K{d) to or to oo. 

Step 2: Propagation of the decomposition into solitons 

Here, we are going to use essential facts of the theory of "solitons", namely that under 
the flow of equation (2.3) and uniformly with respect to r close to ro, this decomposition 
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is stable in time as time increases and that there arc no collisions between solitons. More 
precisely, the idea is to use the equation (2.3) satisfied by Wr to propagate this decompo- 
sition from s = Lm+i to | log |r — ro|| + L where L is large, and prove (roughly speaking) 
the following (see Proposition 3.1 in [18] for a precise statement): 



sup 

Lm+i<s<\ log \r-ro\\+L 



Wris) 
dgWris) 







(4.26) 



as Lm+i — )• oo, L — )• oo and r — )• ro for some parameters {di{s), t'i(s)). 
Let us remark that we can reduce to the case 

ro = r(ro) = 0, 

provided that we change equation (1.4) by the following: 

—drU + \u\'^~^u for all r > — ro- 



o2 r,2 ^ (^ 



This step involves two techniques: 

- a modulation technique of the solution around the sum of solitons. The one-dimensional 

case is treated in Section 2 in [18] which happens to be independent of the equation, hence 
it holds in the radial case with no modifications. More precisely. Section 2 depends only 
on the solitons (4.25) which are the same in one space dimension and in higher dimensions 
in the radial case. 

- the proof of the stability of the decomposition into a decoupled sum of solitons, performed 
in Section 3 of [18] in the one-dimensional case. If many estimates are independent from 
the equation, it is natural that some parts slightly change, because they use the equation 
in similarity variables. It is important to note that even though equation (2.3) depends on 
the considered point r and on time s, the difference with the one-dimensional case comes 
only from one term whose effect is uniformly bounded thanks to (3.8). The estimate (3.8) 
makes it easy to control the effect of the additional term in the adaptation of Appendix 
C of [18], where we project the linearization of equation (2.3) around the sum of solitons. 
More precisely, as we did in Proposition 3.1, we have to add the term Ce~^ to the right- 
hand side of the four differential inequalities of Lemma C.2 in [18]. Accordingly, we have 
to mention that the statement of Claim 3.8 in [18] slightly changes by adding the term 
e~^/S* to the right-hand side of the differential inequalities involving hi and /i2, and also 
by mentioning that 5*{Lm) G (0, 1) and not just that (5*(L^) > 0. 

Step 3: Trapping near one soliton and conclusion 

Following (4.26), it happens that at time s = |log|r — ro|| -|- L, all the solitons 
K*{di{s)^i>i{s)) for i = 2,...,k become small (or vanish) for L large and |r — roj small 
(see Claim 3.4 in [18] for a precise statement), so only the first soliton is left in (4.26). 
Since 



Vr G [- 



and s > ss{ro), E{wr{s)) > E{kq) - Cz{rQ)e 



from (i) of Corollary 7, we see that the first soliton has to be a pure soliton of the form 
— Ac((i^,0) given in (3.1), for some explicit d\ = d\{r), leading to the following estimate: 



dsWr{s*) 



+ 



K{dl) 





< eo where s* = \ log |r — rojj -|- L 
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and eo = eo(^) > is defined in the trapping result stated in (ii) of Theorem 2. Here, there 
is a dehcate point in the adaptation, since the statement of Theorem 2 is different from 
the one-dimensional case, in the sense that we need to apply it uniformly for r > po = 
which is the case whenever \r — ro| is small. Applying that trapping result, we derive two 

facts: 

- the point r is non characteristic, hence ro is an isolated characteristic point (this is the 
conclusion of (i) in Theorem 8); 

- the slope T'{r) satisfies | argthT'(r) — argthJ^(r)| < Ceo, which gives by integration the 
desired estimate in (ii) of Theorem 8. 
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